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ABSTRACT 

The abundance of galaxy clusters is in principle a powerful tool to constrain cosmological pa¬ 
rameters, especially Q m and <X8, due to the exponential dependence in the high-mass regime. 
While the best observables are the X-ray temperature and luminosity, the abundance of gal¬ 
axy clusters, however, is conventionally predicted as a function of mass. Hence, the intrinsic 
scatter and the uncertainties in the scaling relations between mass and either temperature or 
luminosity lower the reliability of galaxy clusters to constrain cosmological parameters. In 
this article, we further refine the X-ray temperature function for galaxy clusters by Angrick 
et al., which is based on the statistics of perturbations in the cosmic gravitational potential 
and proposed to replace the classical mass-based temperature function, by including a re¬ 
fined analytic merger model and compare the theoretical prediction to results from a cosmo¬ 
logical hydrodynamical simulation. Although we find already a good agreement if we com¬ 
pare with a cluster temperature function based on the mass-weighted temperature, including 
a redshift-dependent scaling between mass-based and spectroscopic temperature yields even 
better agreement between theoretical model and numerical results. As a proof of concept, in¬ 
corporating this additional scaling in our model, we constrain the cosmological parameters 
I2 m and <x 8 from an X-ray sample of galaxy clusters and tentatively find agreement with the 
recent cosmic microwave background based results from the Planck mission at lcr-level. 

Key words: cosmology: theory - methods: analytical - dark matter - cosmological parame¬ 
ters - galaxies: clusters: general 


1 INTRODUCTION 

Since the evolution of galaxy clusters is mainly driven by gravity 
and dark-matter (DM) physics, galaxy clusters provide in principle 
reliable information about cosmological parameters, especially on 
the density of matter in the Universe and the amplitude of initial 
density perturbations, expressed by fl m and erg, respectively. The 
usual procedure is to fit a theoretically motivated X-ray temperature 
function to X-ray samples obtained by analysing surveys conducted 
with X-ray telescopes like XMM Newton and Chandra. 

The statistics of DM haloes, however, could conventionally 
only be analytically predicted based on masses, but measurements 
yield usually the X-ray temperature and luminosity. This implies 
that scaling relations between mass and other observables have to 
be calibrated relying on other cluster observables like e.g. weak- 
lensing shear, which reduces the power of the theoretical prediction 
by introducing additional sources of scatter. Additionally, the mass 
of a galaxy cluster is a global quantity and hence ill-defined since 
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galaxy clusters have neither a fixed boundary with a regular shape 
nor a well-defined centre. 

Angrick & Bartelmann (2009) developed an approach to de¬ 
rive the X-ray temperature function from the statistics of perturba¬ 
tions in the cosmic gravitational potential and the relation between 
potential depth and X-ray temperature. Using an extension of the 
ellipsoidal-collapse model of Bond & Myers (1996) by Angrick & 
Bartelmann (2010), Angrick & Bartelmann (2012) refined their X- 
ray temperature function by replacing spherical- with ellipsoidal- 
collapse dynamics and including an analytical merger model that 
takes the temporal temperature increase due to mergers statistically 
into account. 

The goal of this paper is to further refine the merger model 
and analyse how the agreement between their analytical prediction 
and numerical realizations of the temperature function depends on 
the temperature definition used when setting up the halo catalogue 
inferred from a numerical simulation. We also want to analyse how 
constraints on f2 m and erg using the potential-based temperature 
function are affected if differences between the theoretically moti¬ 
vated X-ray temperature that enters via the virial theorem and tem- 
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peratures that are actually inferred from measurements with X-ray 
telescopes are taken into account. 

The structure of this paper is as follows: In Sect. 2, we summa¬ 
rize the main ingredients of the cluster temperature function built 
upon the statistics of the cosmic gravitational potential including 
the ellipsoidal-collapse model by Bond & Myers (1996), slightly 
extended by Angrick & Bartelmann (2010), and the refined ana¬ 
lytic merger model based on the model by Angrick & Bartelmann 
( 2012 ). 

In Sect. 3, we give an overview of the fully hydrodynamical 
numerical simulation we base our analysis upon and introduce the 
temperature definitions we work with in the remainder of this work. 

We contrast our theoretical model to results from this numer¬ 
ical simulation for different temperature definitions in Sect. 4 and 
analyse how taking into account a relation between these differ¬ 
ently defined temperatures improves the agreement between semi- 
analytic model and simulations. 

In Sect. 5, we discuss how we can use these findings to con¬ 
strain the cosmological parameters fl ra and cr 8 from a fit to an X- 
ray sample by Vikhlinin et al. (2009) using the C statistic by Cash 
(1979). 

We present our constraints on both parameters in Sect. 6 and 
quantify how a relation between the X-ray temperature from the 
model and the one inferred from measurements affects constraints 
on both parameters. Additionally, as a proof of concept, we com¬ 
pare to recent results from the Planck mission. 

Finally, we give a short summary in Sect. 7 and provide an 
outlook on how to proceed further with our potential-based tem¬ 
perature function. 


2 THEORETICAL BACKGROUND 

In this section, we will give short overviews of the main ingredients 
that are needed to construct the X-ray temperature function without 
reference to mass. 


2.1 The ellipsoidal-collapse model 

Here, we will present the essential steps towards the ellipsoidal- 
collapse model used later in this work. For detailed calculations, 
refer to Bond & Myers (1996) and Angrick & Bartelmann (2010). 

Let a,- = /?;/R p k with 1 < i < 3 be the dimensionless principal 
axes of a homogeneous ellipsoid, where the Rj are the ellipsoid’s 
physical semi-major axes, and R p ^ is the size of a spherical top-hat 
corresponding to a mass M = (47t/3)pb-R 3 k with the cosmological 
background density pb. 

Their evolution is described by the following three coupled 
differential equations, 


d 2 a,- 

da 2 


1 E'(a) 

a E(a) 


da, 

da 


3a„ 


2a 5 E 2 (a) 


C,(a) ■ 




a 2 E 2 (a) 


oi = 0 , ( 1 ) 


where a is the scale factor of the Universe, E(a) its expansion func¬ 
tion, and E'(a) its derivative with respect to a, f2 m and fl A are 
the dimensionless density parameters of matter and the cosmologi¬ 
cal constant, respectively, in units of the critical density today and 
C, = (1 + 8)/3 + bjl + T ext ,„ 

Here, S = a 3 / ( 0 ^ 303 ) - 1 is the density contrast of the evolv¬ 
ing ellipsoid, b t is the ith component of the internal shear given 


by 


bj(a) = ai{a)a 2 (a) 03 (a) 


I 


dr 


0 [a 2 (a) + 1] UlM(a) + 1] 1/2 3 


( 2 ) 

and the Text,; denote the components of the external shear, which 
we approximate by the so-called hybrid model. 


Text,/(A) — 


^b,(a) 

D+(a ) 

D+(a ta ,,) 


if a < a ta ,, , 

^ext,i(^ta,z) else, 


(3) 


where a ta ,, is the scale factor of turn-around of the ith axis, D + 
is the linear growth factor of matter perturbations, and the T; are 
the eigenvalues of the Zel’dovich deformation tensor (Zel’dovich 
1970). Using the hybrid model, we take into account that the ellip¬ 
soid’s evolution is initially tightly bound to its vicinity until each 
axis finally decouples at its turn-around. 

The initial conditions of equation (1) are given by the 
Zel 'dovich approximation. 


a i( a o) = flo[l - Mao)l 


da, 

da 


ao 


1 - T,(a 0 ) - 


dln£) + 
din a 


Mao) ~ 

a 0 


1 - 2T,(a 0 ), 


(4) 

(5) 


since D + (ao) ~ ao for a small initial scale factor a 0 chosen to be 
a 0 = 2 x 1CT 3 . Note that in this ellipsoidal-collapse model, the 
reference frames of both the ellipsoid and the gravitational shear 
coincide. 

The initial values of the T, at a 0 are given by 


_ 6 a(M) 

3 VTOtt’ 



S o-(M) 

3 yion’ 


(6) 


where a(M) is the square root of the matter power spectrum’s 
variance filtered with a circular top-hat function W M on the scale 
R(M) = [2GM/(fi m tf 2 )] 1/3 given by 


^(M) = f g ^k 2 Ps(k)W 2 M (k). (7) 

Here, M is the halo’s virial mass, G is the gravitational constant, 
Ho is the Hubble constant, Pg(k) is the matter power spectrum and 
W M (k) is the Fourier transform of W M - 

To derive equation (6), we implicitly used that the most prob¬ 
able values for the initial ellipticity and the prolaticity are the best 
choices to describe a statistical average of haloes with mass M. 

According to equation (1), the evolution of the three axes a, 
continues until the smallest axis a; finally collapses first. At that 
point, the halo’s density would formally be infinitely large, and the 
evolution of the other two axes could no longer be followed. But 
physically, the ellipsoid’s collapse should be stopped before due to 
virialization. 

A proper virialization condition can be derived for each axis 
from the tensor virial theorem. It reads 



1 / 3f2„ 


a 2 E 2 (a) \ 2a 3 


C, — 12/ 


and a\ < 0. 


(8) 


If the former expression is fulfilled for the ith axis, its collapse is 
stopped by hand, i.e. a,- is set constant and a' to zero for the further 
evolution. The ellipsoid is considered to be completely collapsed 
when the condition (8) is fulfilled for the largest axis a 3 . This de¬ 
fines the virialization scale factor a v . 

The critical overdensity 8 C and the virial overdensity A v with 
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respect to the critical density are given by 
D+(a v ) 3 


<5 C = —— y x and a v = 

D + (a 0 )^ 


a v D m ((i v ) 


ai(a Y ) a 2 (a v ) a 2 (a v ) 


(9) 


respectively. Note that both 6 C and A v are mass dependent in this 
model unlike results from the spherical-collapse model! 


2.2 The X-ray temperature function without reference to 
mass 


In the following, we will shortly describe how to evaluate the clus¬ 
ter temperature function derived from the statistics of minima in 
a Gaussian random field and incorporating ellipsoidal-collapse dy¬ 
namics. For a detailed derivation, see Angrick & Bartelmann (2009, 
2012 ). 

Based on the statistics of Gaussian random fields (Bardeen 
et al. 1986), the number density of minima in the cosmic gravita¬ 
tional potential as a function of the linear potential depth ®| and the 
potential’s Laplacian AO can be derived analytically and is given by 


n(®i, AO) = 


2407t 3 CTj -ft5y 


{Fi + F 2 )ex p 


(2cr 2 A® + cr 2 ®i) O, 


2y 


with 


F[ = 2tr 2 (5AO 2 - lbcrj) exp 
+ cr 2 (l55A® 2 + 32cr 2 ) exp 


(fiercer 2 - 5crj) AO 2 


2<r 2 y 

(9cr 2 cr 2 - 5o-j) AO' 


8o-;y 

cr 2 A0 2 


F 2 = 5 VTOtiAO (AO 2 - 3cr 2 ) exp 

erff^Uerff^ 

V 2 V2cr 2 ) \flcr 2 ) 


( 10 ) 


( 11 ) 


( 12 ) 


Here, the <Tj with 0 < j < 2 are the spectral moments of the poten¬ 
tial power spectrum Pq,(k) defined as 


cr i = 


= f ^k 2+2i P tb (k)Wl Arb (k), (13) 


and y = cr 0 cr 2 - cr\. The window function in the former equation is 
given by 


VVcd, ao(^c) 


5 [3 sin u - u (3 + m 2 ) cos «] 
2 u 5 


(14) 


where u = kR with R = V _ 2®i/A®, which is the Fourier transform 
of the functional form of a homogeneous sphere’s gravitational po¬ 
tential. 

The lower integration boundary £ min in equation (13) is chosen 
for each pair (®i, AO) such that the number density (10) is maxi¬ 
mized. This effective high-pass filter removes large modes of the 
potential and also of its gradient so that the potential of an object is 
defined with respect to its local environment and small objects with 
VO t 0 are brought to rest and are therefore counted correctly. 
These objects correspond to minima in the non-linearly evolved 
potential which, however, are not present in the linearly evolved 
one. 

The potential power spectrum P®(k) is related to the matter 
power spectrum Ps(k) by 

9 Cl 2 Hi 

p *(k)= --^-}Pdk). (15) 

4 a kr 


To obtain the number density of minima in the gravitational poten¬ 
tial that belong to collapsed structures, we have to integrate equa¬ 
tion (10) over AO accordingly, 


«(®i) = 



dAOn(®,, A®) 0 h [A® - A® c (®i, A®)], 


(16) 


where 0 H is Heaviside’s step function and the critical Laplacian 
A® c is given by 


A® c (®i,A®)= 


<5 c (®i. A®) 


(17) 


Here, 6 C is dependent on both ®i and A® through equations (6) and 
(7) by setting again R = V~2®i/A® instead of R(M). 

Since galaxy clusters are highly non-linear objects, we have 
to relate the linear potential, ®i, to a non- linear one, ®ni. This is 
realized by comparing the potential in the centre of a homogeneous 
ellipsoid at the time of virialization to the linearly propagated po¬ 
tential. The result is 


®ni _ _Ov_ D+(a 0 ) r _dr_ 

®i -do £>+(«v) 0 ^a 2 (a v ) + r] |a;(a v ) + r] [a^(a v ) + r] 

(18) 

where <5 0 is the initial overdensity inside the collapsing halo chosen 
such that the last axis virializes at a v . Note that the ratio ® n i/®i de¬ 
pends on both ®i and A® through the axes a, since the ellipsoidal- 
collapse model is initialized with the scale R = sj- 2®i/A®. How¬ 
ever, since we need to relate a single linear potential to a non-linear 
one, we marginalize over the dependence on A® by 


(®l)A#(®nl) - 


ff dA® ®i n(® 1 , A®) 0 h [A® - A® c (®i, A®)] 
dA® n(®i, A®) 0 h [A® - A® c (®!, A®)] 

(19) 


where ®j is a function of ® n [ and A® via equation (18). 

The non-linear potential depth can be related to an X-ray tem¬ 
perature T via the virial theorem if we average over sufficiently 
many particle orbits. For particles near the centre, this yields 


k^T = -y7u p ®„i, (20) 

where k B is Boltzmann’s constant, m p is the proton mass and p 
represents the mean molecular weight and is defined as 

P = ( 2X P + 1 7 p) • (2D 


where X v = 0.76 and F p = 1 - X p are the primordial hydrogen and 
helium mass fractions, respectively, so that p x 0.59. 

Summarising the above relations, the number density of gal¬ 
axy clusters as a function of their X-ray temperature is given by 


( 20 ) ( 19 ) 

n Va {T) AT = n(T J ® nl -> <®,) A *) 


( 16 ) 


d(®l>AO d® n | 

d®„, "dF 


dT, 


( 22 ) 


where the derivative d(®i) A< i ) /d® nl has to be calculated numerically 
from equations (18) and (19). 


2.3 The influence of mergers on the temperature function 

So far, we have only taken into account virialized structures for 
which the relation (20) holds. However, due to hierarchical struc¬ 
ture formation, smaller DM haloes merge to form larger haloes. 
Since these mergers induce a rise in the X-ray temperature func¬ 
tion (see e.g. Randall, Sarazin & Ricker 2002), we will include 
this effect statistically in our framework by a simple parameter-free 
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merger model. The following model slightly extends the model by 
Angrick & Bartelmann (2012) to account for the change of the clus¬ 
ter temperature function with redshift. 

Starting from the number density of virialized galaxy clusters 
at a given temperature and redshift, n yk (T, z), we calculate two cor¬ 
rection terms. First, we add the clusters that reach a temperature T 
only due to a temperature boost AT induced by a merger, and sec¬ 
ond, we subtract those that would have a temperature T if they were 
virialized, but have a temperature higher than T due to a merger. 

Denoting the first contribution by n + (T) and the second by 
n_(r), the final halo population n(T,z ) can be modelled as 

n(T, z) = n vir (r, z) + n + (T, z) - n.(T, z), (23) 


where n + (T, z) is given by 


r»Z+tAz s * 00 

n+(T,z)= | dz' I dA M I dMn[7 vir (M,z'),z'] 

Jz Jo Jo 

x p(M, AM, z ) & d [T - T vir (M, z) - AT(M, AM, z')] 
x 0 H (M - AM) (24) 


and n_(T, z) by 

f Z+A z r'M(T,z') 

dz' I dAMn(T,z')p[M(T,z'),AM,z'l 

(25) 

In the two former equations, there are various quantities that have 
to be defined in the following. 

The merger rate p(M,AM,z)dAMdz yields the number of 
mergers of haloes with mass M with other haloes in the mass range 
[AM, AM + dAM] and in the redshift range [z, z + dz] and is given 
by 


p(M, AM, z) 


1 

V27t 


Si 

Si(Si-S 2 ) 


3/2 

exp 


dS 2 dio 
dAM dz 


ur(Si-S 2 ) 

2SiS 2 


(26) 


(Lacey & Cole 1993), where S i = a 2 (M), S 2 = cr(M + AM) (cf. 
equation 7), and ai = 8 c (z)/D + (z ). Here, 8 c (z) is the critical over¬ 
density of the spherical-collapse model. Heaviside’s 0-function in 
equation (24) ensures that M > AM so that mergers between two 
objects with masses M and AM are only counted once by demand¬ 
ing that the mass ratio between main and merging object is always 
larger than unity. 

The redshift interval Az is connected to the sound-crossing 
time r sc by 


_ R _ 1 f dz' 

tsc ~ c s ~ T/ 0 J E(z')(l + z')’ 


(27) 


where 


R(M) = 


/ 3Maifl 2 fl 3 \ 1/3 

l 47tp b ) 


and c s 



are the cluster’s mean radius and the sound speed, respectively. 

The theoretical temperature-mass relation for virialized ellip¬ 
soidal haloes is given by 


k B T yk (M,z) = pm p 


Q m B 2 G 2 M 2 \ 1/3 


128 


) I 


dr 


0 Yll l [a 2 k (M,z) + T\ 


1 / 2 ’ 


(29) 

Dirac’s delta-distribution 5 D in equation (24) ensures that only 
those clusters contribute to the integral whose temperature T is the 


sum of their temperature based on virial equilibrium, T vil , and a 
temperature increase AT due to the merger given by 


k B 


AT(M, AM, z) = - 


2 pm p AM (®>(M, z) 

3 M 


(30) 


where (®)(M, z) is the mean potential on the surface of an ellipsoid 
with mass M, virialized at redshift z, 


<0>(M,z) = - 


3 MG 


f 


dr 


4R pk (M) (J 0 nL, [al(M,z) + T\ 


1/2 


\a 2 {M,z) + r] nLi [a 2 k (M,z) + r] 


1/2 I 


(31) 


n_(T) can be calculated accordingly from equation (25) using equa¬ 
tion (26) for p(M, AM, z), equation (27) for Az, and inverting equa¬ 
tion (29) to arrive at M(T, z). Note that in equation (25), T is simply 
the cluster’s virial temperature T yk . 


3 THE SIMULATION 

Hydrodynamical simulations are important numerical tools to de¬ 
scribe in detail the evolution of cosmic structures in the Universe 
taking into account gas physics. In this section we describe in de¬ 
tail the numerical set-up adopted in this work. In particular, we will 
describe the simulations used, how haloes were extracted and how 
we built the halo catalogues with their X-ray properties. 

The hydrodynamical simulation used in this work is part of 
a bigger set of simulations (BINGS, Baryons in non-Gaussian 
.simulations) originally proposed to study the effect of baryons in 
cosmologies with non-Gaussian initial conditions. They represent 
an extension of the DM only simulations analysed to study sev¬ 
eral aspects of the large-scale structures of this kind of cosmologies 
(Grossi et al. 2009). In addition, the DM-only simulations with the 
same cosmological parameters and box size have been already em¬ 
ployed for studies on the SZ effect and X-ray emission (Roncarelli 
et al. 2010). For this work, we use exclusively the simulation run 
with Gaussian initial conditions, also in the light of recent observa¬ 
tions of the cosmic microwave background (CMB) by the Planck 
mission (Planck Collaboration XIII 2015). 

The simulations are based on the parallel cosmologi¬ 
cal TreePM-SPH code Gadget-2 (Springel 2005) using an 
entropy-conserving formulation of smoothed-particle hydrody¬ 
namics (SPH; Springel & Hernquist 2002). The simulation run in¬ 
cludes radiative cooling, star formation with associated feedback 
processes (Springel & Hernquist 2003) and heating by a uniform, 
time-dependent ultraviolet background. The simulation also takes 
into account a multiphase model for star-formation and feedback 
processes due to supernovae-driven galactic winds (Springel & 
Hernquist 2003; Di Matteo et al. 2008). 

The numerical set-up assumes a box of 1.2 h~ l Gpc comoving 
with 2 x 960 3 DM and gas particles. The underlying cosmological 
model is a standard A cold dark matter model with the cosmologi¬ 
cal parameters derived by the 5 yr WMAP results. We report them 
in Table 1. DM particles have a mass of 1.171 x 10 11 IP 1 M G , and 
gas particles have a mass of 2.38x 10 10 h~' M e (at the redshift of the 
initial conditions, z; n i = 60). The gravitational force is computed as¬ 
suming a Plummer-equivalent softening length of e = 25 IP 1 kpc. 
The evolution of particles is followed till the present time (z = 0). 
To create the initial conditions of our simulation, the transfer func¬ 
tion by Eisenstein & Hu (1998) was used. 
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2-0 2 = 0.507308 2= 1.00131 





2= 1.50557 



7™ (keV) 


2 = 2.04689 



Figure 1. The analytic temperature function n(T) including merger effects (equation 23) and the temperature function inferred from the numerical simulation 
based on 7',“™ at five different redshifts. 


Table 1. Set of cosmological parameters at z = 0 employed in the simu¬ 
lations. In the lower part additional characteristics of the simulations are 
reported. 


Parameter 

Symbol 

Value 

Hubble parameter (100 km s _1 Mpc -1 ) 

h 

0.72 

Amplitude of fluctuations at 8 /z _1 Mpc 

0-8 

0.8 

B ary on density 

n b 

0.044 

Total matter density 

Dm 

0.26 

Dark energy density 

Da 

0.74 

Initial redshift 

£ini 

60 

Cubic box length (/i -1 Mpc) 

A)OX 

1200 

Total number of particles 

^part 

2 x 960 3 

Mass of DM particles (h~ l M 0 ) 


1.171 x 10 11 

Mass of gas particles at zj n j (h~ l M©) 

™ g as 

2.38 x 10 10 


Haloes were initially identified with a Friend-of-Friend algo¬ 
rithm (Davis et al. 1985) with linking length b = 0.2 times the mean 
interparticle distance and in a second step, using the SUBFIND 
algorithm (Springel et al. 2001; Dolag et al. 2009), we evaluated 
spherical-overdensity masses for each halo centred at the deepest 
potential point. We computed two different overdensities, A = 200 
and A = 500 times the critical density. The corresponding radii are 
/?200 and R;oo, respectively. We extracted our catalogues at differ¬ 
ent redshifts, namely z = 0, z ~ 0.5, z ~ 1, z ~ 1.5, and z ~ 2, 
and considered only objects with at least 100 DM particles to in¬ 
corporate only haloes in the further analysis whose X-ray profiles 
are resolved well enough.. 

To evaluate the X-ray properties of each halo, we took into 
account all its gas particles, computed the relevant quantities for 
each individual particle and summed them up to have the integrated 
value for the object considered. For what concerns us, we will limit 
our discussion to the evaluation of the halo temperature. Starting 


from the internal energy u of a given SPH particle, we evaluate its 
temperature with the following relation, 

2 

Tj = —gm p i(j. (32) 

3 k B 

For a given object we can define three different temperatures: 
the mass-weighted one, T mw , (Kang et al. 1994; Bartelmann & 
Steinmetz 1996; Mathiesen & Evrard 2001), the emission-weighted 
one, T ew , (Kang et al. 1994; Mathiesen & Evrard 2001), and the 
spectroscopic-like temperature T s i (Mazzotta et al. 2004). They are 
defined as 


f d VpT 
/Wp ’ 
f d Vp-T 
f dVp 2 ’ 
f d Vp 2 T‘ /4 
f dVp 2 T- 3 ' 4 ' 


(33) 

(34) 

(35) 


respectively, where p is the density and V is the integration volume 
(in this case the volume of the halo). In the following, we will con¬ 
sider a spherical integration volume with either R 2 oo or Rsoo as a 
radius and indicate the corresponding temperature as, e.g., T 2< f or 
T™, respectively. 

Since the halo density is not simulated as a continuous quan¬ 
tity, we can replace the integration with a sum over all particles 
such that the operative definitions of the three different tempera¬ 
tures become 


T* 



(36) 



Z, WiPiTi 

(37) 

Zi m iPi 


Z, 1’iiPiTf 4 

(38) 

Z,m,p,rr 3/4 ’ 
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z = 0 



Z = 0.507308 



z — 1.00131 



Z = 1.50557 



Z = 2.04689 



Figure 2. Same as Fig. 1 , but the temperature function inferred from the simulation is based on T ^ 00 ' cut 


where m, and p, are the mass of the fth particle and the SPH estimate 
of the density at its position, respectively. 1 

The mass-weighted temperature (equation 33) has a relevant 
physical meaning since it is proportional to the total thermal en¬ 
ergy of the cluster. However, it may differ significantly from what 
an observer would measure via an X-ray spectral analysis since 
the bremsstrahlung emissivity is proportional to the square of the 
gas density: the emission-weighted temperature (equation 34) takes 
this bias into account. The spectroscopic-like temperature (equa¬ 
tion 35), introduced by Mazzotta et al. (2004), besides the emis¬ 
sivity bias, accounts also for the gas temperature gradient along the 
line of sight together with the shape of the energy response function 
of Chandra and XMM-Newton X-ray detectors. 2 


4 ANALYTICAL MODEL VERSUS SIMULATION 

In the following, we want to compare the results of our analytical 
model from Sect. 2 with results from the simulation presented in 
the previous section based on the two temperature definitions 7j 2 ®’ 
and 'T’™ 1 ' cut in the analysis, where ‘cut’ in the superscript refers to 
the fact that the volume corresponding to the inner 15 per cent of 
the sphere’s radius is cut out from the integration to obtain T 300 . 
The reason for cutting out the inner part will become clear shortly. 

In Fig. 1, we show the X-ray temperature function based on 
f™ together with our semi-analytical prediction for five different 
redshifts. The theoretical curve matches the numerical results from 

1 Note that there are more complicated definitions of the emission- 
weighted temperature, involving in its definition the cooling function 
(Bryan & Norman 1998; Frenk et al. 1999; Muanwong et al. 2001; Borgani 
et al. 2004). However, since we will not use the emission-weighted tem¬ 
perature in the further discussion, we limit ourselves to the more simplified 
definition (34) here. 

2 When calculating T ew and 7(], it is necessary to remove cold SPH parti¬ 
cles (7) < 0.5 keV) from the computation to exclude the contribution of the 
cold dense phase of the ICM that has negligible emissivity. 


the simulation very well for lower temperatures. At temperatures 
>4-5 keV, however, the theoretical prediction lies above the tem¬ 
perature function from the simulation. A reason for the latter might 
be that haloes with a high X-ray temperature and therefore also a 
large mass are not well resolved in the simulation due to the limited 
number of large-scale modes in the simulation. 

Since the X-ray temperature in our model is derived from the 
thermal energy of the virialized cluster and the mass-weighted tem¬ 
perature is proportional to the cluster’s thermal energy, 7’ 2 ® 1 is the 
definition used in the simulation that matches best the one in our 
model and hence is the appropriate quantity in this case. 

In Fig. 2, we show the X-ray temperature function based on 
y500, cut as j n f erre( j f rom the simulation together with our semi- 
analytic prediction since X-ray observers often fit the spectrum of 
a galaxy cluster in the interval [0.15, l]/?soo- The theoretical predic¬ 
tion underestimates the number density of haloes at temperatures 
<4-5 keV, whereas for higher temperatures the theoretical curve 
lies above the results from the simulation at low redshifts. This 
might be another hint at missing high-mass objects in the simula¬ 
tion and incomplete statistics. Overall, the theoretical curve seems 
too flat compared to the numerical one inferred from the simula¬ 
tion. 

To relate the temperature that is measured by observers (~ 
Tf °- cut ) to the one that our model is based on (~ T 2 “), we fitted a 
relation for each of the five different redshifts of the form 

y200 \ b 

3”keV/ ■ (39) 

based on the cluster catalogue extracted from our simulation. The 
fitting parameters T 0 and b therefore vary with redshift. We used 
only those clusters for the fit that fulfil both T > 1 keV and M dm > 
1.17 x 10 13 h~ l M 0 for both r e [0.15,1 ]/?soo and r < R 200 , where 
M dm is the mass of the cluster’s DM component. The cut in mass 
corresponds to the requirement of having a halo that consists of at 
least 100 particles. The fitted functions (39) together with the data 
points the fits are based on are shown for various redshifts in Fig. 3. 

The resulting fitting parameters To and b can be found in Ta- 
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Figure 3. Relation between T™ and T 3011 ' cut inferred from the halo catalogue including fits of the form (39) for five different redshifts. 


Table 2. The best-fitting parameters To and b as well as their errors A To 
and A b, respectively, as a function of redshift z. 


z 

To (keV) 

b 

AT 0 (keV) 

A b 

0 

3.209 

0.9012 

0.003 

0.0013 

0.507308 

3.085 

0.9114 

0.005 

0.0020 

1.00131 

3.061 

0.9399 

0.007 

0.0029 

1.50557 

3.051 

0.9565 

0.012 

0.0048 

2.04689 

3.075 

0.9717 

0.023 

0.0087 


ble 2. Except for z = 0, the parameter T 0 is very close to 3 keV. 
Additionally, the higher the redshift, the closer to unity is the pa¬ 
rameter b. The rms of the deviations of individual clusters from the 
fitting relation (39) is with 0.21 keV largest for z = 0 and with 
0.12 keV smallest for z = 2 and hence relatively small. 

Based on the relation (39), we recalculated our theoretical pre¬ 
diction by scaling the temperature function accordingly, i. e. given 
7^0°, cut , we calculated T™ and used this temperature as input for 
equation (22). In doing so, an additional factor |dr™/dj r s 5 1 00 ' cut | had 
to be taken into account. The comparison between the rescaled the¬ 
oretical temperature prediction and the one inferred from the simu¬ 
lation based on T 300- ™' can be found in Fig. 4. 

After a temperature rescaling in the theoretical prediction for 
the X-ray temperature function, the agreement between the theo¬ 
retical predictions and the numerical results based on pf°- cut is 
much better for all redshifts. Thus, by taking the redshift-dependent 
temperature conversion (equation 39) into account, our theoreti¬ 
cal model based on the statistics of gravitational-potential pertur¬ 
bations is consistent with numerical results that mimic as well as 
possible the X-ray temperature function inferred by real observa¬ 
tions. 


5 CONSTRAINING COSMOLOGICAL PARAMETERS 

Based on the results from the previous section, we want to deter¬ 
mine the cosmological parameters <r 8 and fl m from a cluster sample 
by Vikhlinin et al. (2009), which was also considered in the Planck 
2013 analysis of the Sunyaev-ZeTdovich (SZ) effect (Planck Col¬ 
laboration XX 2014), using our theoretical model for the X-ray 
temperature function together with a proper temperature scaling de¬ 
pendent on redshift based on equation (39). The statistical analysis 
is very similar to the one already presented in Angrick & Bartel- 
mann (2012). 

This section is meant as a proof of concept. To reliably con¬ 
strain cosmological parameters in the future, one should stick to 
empirical relations between the potential of a cluster and its tem¬ 
perature, e.g. by combining data from gravitational lensing and X- 
ray temperature measurements, and not rely on relations inferred 
only from hydrodynamical simulations. See also the end of Sect. 7 
for more details. 


5.1 The sample 

The sample by Vikhlinin et al. (2009) consists of two subsamples, 
one at high and one at low redshift, based on ROSAT PSPC All- 
Sky (RASS) and 400 deg 2 data. The low-redshift sample is based 
on archival Chandra data and consists of 49 clusters with flux / > 
1.3xl0 -n ergs -1 cm -2 in the 0.5-2 keV band from several samples 
of RASS with a total area of 8.14 sr (26722 deg 2 ). The redshift 
coverage is 0.025 < z < 0.25 with (z) « 0.05, and temperatures 
are in the range 2.61 keV < T < 14.72 keV as inferred from their 
spectra measured in multiple annuli with Chandra. 

The high-redshift sample consists of 36 clusters from the 
ROSAT 400 deg 2 survey (Burenin et al. 2007) further analysed with 
Chandra in the redshift range 0.35 < z < 0.9 with (z) ~ 0.5 
and a redshift-dependent flux limit in the 0.5-2 keV band. For 
z > 0.473, the limiting flux is 1.4 x 10 -13 erg s -1 cm -2 , whereas for 
z < 0.473, the flux limit corresponds to a minimal X-ray luminosity 
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z = 0 z = 0.507308 z = 1.00131 





Z = 1.50557 z = 2.04689 




Figure 4. Same as Fig. 2, but the theoretical temperature function was rescaled according to the redshift-dependent relation between T™ and 7^°°’ cut inferred 
from Fig. 3. 


of L x , min = 4.8x 10 43 (1 + z) L8 erg s _l . The temperatures of the clus¬ 
ters are in the range 2.13 keV < T < 11.08 keV as inferred from 
their spectra measured with Chandra in the region [0.15, l]i? 5 oo, 
where R 500 is the radius that encloses a mean overdensity of 500 
times the critical density of the Universe. 

The effective differential search volume dV/dz as a function of 
mass M and cosmological parameters £2 m , £T A and Hq = 72 km s _1 
Mpc -1 for both subsamples was made available in electronic form 
on a grid by A. Vikhlinin. To convert it to a function of temperature, 
we used the best-fitting values of the mass-temperature relation of 
Vikhlinin et al. (2009), 

M 50 o = Mo(^^) E~\z), (40) 

where M 0 = (3.02 ± 0.11) x 10 14 hr 1 M 0 and a = 1.53 ± 0.08. 


5.2 The fitting procedure 


Since the errors on the cluster number counts are Poissonian, we 
used the C statistic of Cash (1979) for unbinned data to find the 
best-fitting values for fi m and cr 8 , assuming a spatially flat uni¬ 
verse, hence fl A = 1 - fl m . In the next section, we compare with 
the results from the 2015 data release of the Planck Collabora¬ 
tion so that we set the baryon density parameter G b h 2 = 0.02225, 
the index of the primordial power spectrum n s = 0.9645 and 
H 0 = 67.27 km s -1 Mpc -1 as inferred from the TT,TE,EE+lowP 
data (Planck Collaboration XIII 2015). To include baryonic effects, 
we used the transfer function by Eisenstein & Hu (1998) when 
computing the power spectrum Ps(k). 

The C statistic is defined as 


C = 


2 





(41) 


where N is the total number of objects expected from the sample 
assuming a theoretical model, and n, is the theoretically expected 
differential number density of the zth cluster in the sample with 


temperature T, and redshift Zi- The sum extends over all sample 
members. 

Although the cosmological parameters in the further discus¬ 
sion will slightly differ from their numerical values in our cosmo¬ 
logical simulation, we assume that the temperature relation (39) 
does not change significantly so that we could still use the best¬ 
fitting values from Table 2. To take the listed uncertainties into ac¬ 
count, we convolved with a normal distribution of the form 


p(T\T™)dT = 


T-T; 


500, cut 


s/2n 


CT j 


■exp 


(r 2°0)l 2 

V mw 'J 


2 cr- 


d T, (42) 


where 7^ 00, cut (L™) is given by equation (39). Both T 0 and b in 
that relation depend on redshift as follows. Let z„ and z n +i be two 
consecutive redshifts for which the former relation can be deduced 
from the simulation. For a redshift z with z„ < z < z n + 1 , the val¬ 
ues of the parameters T 0 and b as well as their errors are linearly 
interpolated from their values at z„ and z n + 1 (cf. Table 2). To prop¬ 
erly take the errors on the parameters T 0 and b into account, the 
standard deviation was set to 

- C w “<® ,431 

due to Gaussian error propagation. Thus, the expected number of 
objects in each subsample is given by 

Mowihigh = £dz£dT dVl °^\ T , z) (44) 

x f 

where the subscripts ‘low’ and ‘high’ denote the low- and the high- 
redshift subsample, respectively. The integral boundaries depend 
on the subsample and are given in Sect. 5.1 for z and T. The in¬ 
tegration over L™ has to be done over the whole valid range of 

P(7lO- 
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Finally, the expected differential number density of the fth 
cluster is simply given by the convolution 


dl4,w|high 

^/',low|high — j_ \*iiZi) 


d z 


(45) 


l‘ T =l 


dTn(T^)p(T\T 2 Z)qmT), 


with 


q(T,\T) = 


f2na-j 


exp 


(T, - Tf 


2 07 


(46) 


where <x, is the measurement error of the ith cluster. 

To jointly fit both the low- and the high-redshift cluster sam¬ 
ples of Vikhlinin et al. (2009), we had to add the two contributions, 
resulting in 



C = 2 


Mow - I> 


+ Mhigh - V In i 


(/.high 


(47) 


We searched for minima of the C statistic as a function of the two 
cosmological parameters fl m and erg, which enter both via n(Tff) 
and the volume factor dV/dz. 

Cash (1979) showed that one can create confidence intervals 
for the C statistic exactly in the same way as it can be done for a 
X~ fit using properties of the x 1 distribution. Following the work 
by Lampton, Margon & Bowyer (1976), intervals with confidence 
y are implicitly given solving 

y = f dr 2 /Cr), (48) 

Jo 

for t, where / is the density of the x], distribution with p degrees 
of freedom determined by the number of parameters. For 68 per 
cent confidence and p = 2, it follows that t = 2.3, while for 95 per 
cent confidence, t = 5.991. Using the minimum of the C statistic, 
C m i„, we could simply calculate the 68 and 95 per cent confidence 
contours by searching for points in the parameter space for which 
C = C m in + 2.3 and C = C m i n + 5.991, respectively. 


6 RESULTS 

In Fig. 5, we compare the 68 and 95 per cent confidence contours 
for the parameters and cr 8 inferred with our semi-analytic clus¬ 
ter temperature function to the respective confidence contours from 
the Planck 2015 data release (Planck Collaboration XIII 2015) 
based on the TT,TE,EE+lowP data both including the temperature 
relation (39) and excluding it by setting p(T\T™) = S D (T - Tff) 
in equations (44) and (45) for the latter case. 

Taking the temperature conversion from T™’ cut to T™ into 
account when fitting our theoretical X-ray temperature function to 
the X-ray sample by Vikhlinin et al. (2009) shifts the confidence 
contours to smaller values for and to larger values for cr 8 com¬ 
pared to the case where such a temperature conversion based on 
the difference between measured and theoretically motivated tem¬ 
perature is neglected. Additionally, the sizes of both the 68 and the 
95 per cent contour are smaller if the temperature conversion is in¬ 
cluded. 

The Planck constraints are in agreement with our results incor¬ 
porating the redshift-dependent temperature conversion (39) at lcr- 
level, whereas the results without the temperature conversion are in 
agreement with the Planck results only at the 2<x-level. Hence, dis¬ 
criminating between the temperature that is actually inferred from 


Figure 5. 68 and 95 per cent confidence contours for the parameters fi ra 
and erg based on our cluster temperature function as inferred from the two 
subsamples by Vikhlinin et al. (2009) excluding (red contours) and includ¬ 
ing (blue contours) the temperature conversion (39). For comparison we 
also provide the purely CMB-based Planck TT,TE,EE+lowP results (Planck 
Collaboration XIII 2015, black contours). 


an X-ray measurement (similar to r^ 00 ’™') and theoretically mo¬ 
tivated temperature consistent with the virial theorem (similar to 
T^j) seems crucial when our potential-based X-ray temperature 
function is used to constrain cosmological parameters. 


7 SUMMARY & CONCLUSIONS 

In the first part of this article, we further refined the X-ray tem¬ 
perature function for clusters based on the cosmic gravitational 
potential by Angrick & Bartelmann (2009, 2012) by incorporat¬ 
ing the redshift evolution of the temperature function for virialized 
structures (equation 22) into our analytical merger model. We then 
compared our theoretical model to X-ray temperature functions in¬ 
ferred from a fully hydrodynamical simulation based on two differ¬ 
ent temperature definitions at five different redshifts: (1) the mass- 
based temperature inside R 2 oo, T™, (2) the spectroscopic like tem¬ 
perature inside Moo with the inner part <0.15 Moo cut out, T ^ 00 ’ cut . 
The main results can be summarized as follows: 

(i) Our theoretical temperature function is in very good agree¬ 
ment with the numerical one based on 7j™ except for relatively 
high temperatures which is presumably due to resolution effects of 
the simulation in the high-mass regime. 

(ii) Compared to the numerical temperature function based on 
j50Q, cut , Qur theoretical X-ray temperature function underestimates 
the abundance of objects over a large temperature range for all red- 
shifts examined. 

(iii) For each redshift z separately it is possible to find a rela¬ 
tively tight relation of the form T™'™ 1 = r 0 (T™/3 keV)*, where 
both To and b are fitting parameters depending on z. The rms of the 
deviations of individual clusters in the simulation from the above 
relation reaches from 0.12 keV at z = 2 to 0.21 keV at z = 0. 

(iv) Scaling the temperature of our theoretical temperature func¬ 
tion according to the above relation for each redshift and then com¬ 
paring it to the numerical temperature function based on 7’ s 5 | 00 ' cut 
yields very good agreement between both functions. 

The second part of this article was meant as a proof of concept. 
Here, we used the redshift-dependent temperature scaling found in 
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the first part for our theoretical model to fit it to an X-ray sample by 
Vikhlinin et al. (2009) and constrain the cosmological parameters 
Q ra and cr 8 . For comparison we also fitted these parameters with¬ 
out taking the aforementioned temperature conversion into account. 
The main results are the following: 

(v) Incorporating the temperature conversion shifts the confi¬ 
dence contours in the fl m -cr 8 plane towards smaller values of fl m 
and larger values of cr s . Additionally, the size of the contours is 
reduced. 

(vi) We find agreement at lcr-level between the purely CMB- 
based TT,TE,EE+lowP results from the Planck 2015 data release 
(Planck Collaboration XIII 2015) and our theoretical model incor¬ 
porating the temperature conversion, whereas the two only agree at 
2 cr-level with each other if the temperature conversion is neglected. 

Concluding, it seems necessary to establish an empirical relation 
between the potential and the X-ray temperature of a cluster since 
relying solely on a scaling inferred from a hydrodynamical simu¬ 
lation may not allow to constrain cosmological parameters reliably 
enough. Combining data from gravitational lensing and X-ray ob¬ 
servations of both relaxed and merging galaxy clusters could be a 
promising way to find such a relation since gravitational lensing 
probes the gravitational potential of a cluster projected along the 
line-of-sight, whereas one can infer the temperature of its intra¬ 
cluster medium (ICM) from X-ray spectra. 

Once such an empirical relation is established, we will use 
our potential-based temperature function on more recent and big¬ 
ger samples of X-ray clusters to constrain O m and <x 8 from cluster 
cosmology with higher accuracy and reliability. However, finding a 
well-defined physical and analytical model for this relation would 
be even more desirable since it would allow directly the modelling 
of the temperature abundance of galaxy clusters, thus avoiding ad¬ 
ditional sources of scatter. 


ACKNOWLEDGEMENTS 

We thank Klaus Dolag for providing the initial conditions of the 
simulation and Margherita Grossi for the code to run them. Simu¬ 
lations and cluster catalogues were created on the Intel SCIAMA 
High Performance Compute (HPC) cluster which is supported by 
the ICG, SEPNet and the University of Portsmouth. 


REFERENCES 

Angrick C., Bartelmann M., 2009, A&A, 494, 461 
Angrick C., Bartelmann M., 2010, A&A, 518, A38 
Angrick C., Bartelmann M., 2012, A&A, 538, A98 
Bardeen J. M., Bond J. R., Kaiser N., Szalay A. S., 1986, Apl, 
304, 15 

Bartelmann M., Steinmetz M., 1996, MNRAS, 283, 431 
Bond J. R„ Myers S. T„ 1996, ApIS, 103, 1 
Borgani S., Murante G., Springel V., Diaferio A., Dolag K., 
Moscardini L., Tormen G., Tornatore L., Tozzi P, 2004, MN¬ 
RAS, 348, 1078 

Bryan G. L„ Norman M. L„ 1998, ApJ, 495, 80 
Burenin R. A., Vikhlinin A., Hornstrup A., Ebeling H., Quintana 
H„ Mescheryakov A., 2007, ApJS, 172, 561 
Cash W„ 1979, Apl, 228, 939 

Davis M„ Efstathiou G., Frenk C. S., White S. D. M„ 1985, Apl, 
292, 371 


Di Matteo T., Colberg I., Springel V., Hernquist L., Sijacki D., 
2008, Apl, 676, 33 

Dolag K., Borgani S., Murante G., Springel V., 2009, MNRAS, 
399, 497 

Eisenstein D. J., Hu W„ 1998, Apl, 496, 605 
Frenk C. S„ White S. D. M„ Bode P. Bond I. R„ Bryan G. L„ Cen 
R., Couchman H. M. R, Evrard A. E., Gnedin N., et al. 1999, 
ApJ, 525, 554 

Grossi M., Verde L., Carbone C., Dolag K., Branchini E., Iannuzzi 
F„ Matarrese S„ Moscardini L„ 2009, MNRAS, 398, 321 
Kang H., Cen R., Ostriker J. P, Ryu D., 1994, ApJ, 428, 1 
Lacey C., Cole S., 1993, MNRAS, 262, 627 
Lampton M., Margon B.. Bowyer S., 1976, ApJ, 208, 177 
Mathiesen B. F., Evrard A. E., 2001. ApJ, 546, 100 
Mazzotta P, Rasia E., Moscardini L., Tormen G., 2004, MNRAS, 
354, 10 

Muanwong O., Thomas P. A., Kay S. T., Pearce F. R., Couchman 
H. M. P, 2001, ApJ, 552, L27 

Planck Collaboration XIII 2015, preprint (arXiv: 1502.01589) 
Planck Collaboration XX 2014, A&A, 571, A20 
Randall S. W„ Sarazin C. L„ Ricker P. M„ 2002, ApJ, 577, 579 
Roncarelli M., Moscardini L., Branchini E., Dolag K., Grossi M., 
Iannuzzi F., Matarrese S., 2010, MNRAS, 402, 923 
Springel V., 2005, MNRAS, 364, 1105 
Springel V., Hernquist L„ 2002, MNRAS, 333. 649 
Springel V., Hernquist L„ 2003, MNRAS, 339, 289 
Springel V., White S. D. M., Tormen G., Kauffmann G., 2001, 
MNRAS, 328, 726 

Vikhlinin A., Burenin R. A., Ebeling H., Forman W. R., Hornstrup 
A., Jones C., Kravtsov A. V., Murray S. S., Nagai D., Quintana 
H„ Voevodkin A., 2009, ApJ, 692, 1033 
ZeFdovich Y. B„ 1970, A&A, 5, 84 


© 0000 RAS, MNRAS 000, 000-000 


